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Abstract 

The qq spectrum is studied within a chiral constituent quark model. It pro- 
vides with a good fit of the available experimental data from light (vector 
and pseudoscalar) to heavy mesons. The new D states measured at differ- 
ent factories are studied. The ++ light mesons are analyzed as qq pairs or 
tetraquark structures. 

I. INTRODUCTION 

The continuously increasing huge amount of data and its simplicity have converted me- 
son spectroscopy in an ideal system to learn about the properties of QCD. Since Gell-Man 
conjecture, most of the meson experimental data were classified as qq states according to 
SU(N) irreducible representations. Nevertheless a number of interesting issues remains still 
open as for example the understanding of some data recently obtained on the B factories or 
the structure of the scalar mesons. The recently measured D s states are hardly accommo- 
dated in a pure qq description. They present a mass much lower than the naive prediction 
of constituent quark models. Besides, the underlying structure of the scalar mesons is still 
not well established theoretically. In a pure qq scheme the light scalars could be identified 
with the isoscalars / (600) and / (980), the isodoublet k(900), and the isovector a (980), 
constituting a SU(3) flavor nonet. However such identification immediately faces difficulties 
to explain for example (i) the / (980) and a (980) mass degeneracy, (ii) why a (980) and 
/o(600) have so a different mass?, (iii) the similar branching ratios of the J/ip — > /o(98O)0 
and J/ip — > /o(980)c<j decays which clearly indicates the existence of strange and nonstrange 
content in the / (980). 

The theoretical tools to determine the properties of mesons are based to a large extent 
on phenomenological models. The study of charmonium and bottomonium made clear that 
heavy-quark systems are properly described by nonrelativistic potential models reflecting 
the dynamics expected from QCD [1]. For heavy quarks to leading order in v 2 it has been 
demonstrated that the interaction can be derived from the theory [2]. The light meson 
sector has been studied by means of constituent quark models, where quarks are dressed 
with a phenomenological mass and bound in a nonrelativistic potential, usually a harmonic 
oscillator [3]. Quite surprisingly a large number of properties of hadrons could be reproduced 
in this way [4]. In this talk we present the meson spectra obtained by means of a chiral 
constituent quark model in a trial to interpret some of the still unclear experimental data 
in the light scalar and D meson sectors. 
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II. SU(3) CHIRAL CONSTITUENT QUARK MODEL 



Since the origin of the quark model hadrons have been considered to be built by con- 
stituent (massive) quarks. Nowadays it is widely recognized that the constituent quark mass 
appears because of the spontaneous breaking of the original SU(3) L <g> SU(3) R chiral sym- 
metry at some momentum scale. In this domain a simple Lagrangian invariant under the 
chiral transformation can be derived as [5] 

L = ^/j(ry M <9 M — MU l5 )ip (1) 

where U 15 = exp(i7r a A a 7s// 7r ). n a denotes the pseudoscalar fields (n,Ki,r]s) with i=l,...,4, 
and M is the constituent quark mass. An expression of the constituent quark mass can be 
obtained from the theory, but it also can be parametrized as M(q 2 ) = m q F(q 2 ) with 



f(q 2 ) 



A 2 



A 2 + q 2 



(2) 



where A determines the scale at which chiral symmetry is broken. Once a constituent quark 
mass is generated such particles have to interact through Goldstone modes. Whereas the 
Lagrangian ip^ij^d^ — M)ip is not invariant under chiral rotations, the one of Eq. (1) is 
invariant since the rotation of the quark fields can be compensated renaming the bosons 
fields. C/ 7s can be expanded in terms of boson fields as, 

C/ 7S = 1 + ^ 5 X a 7r a - -^ttV + - (3) 

Jit ^Jtt 

The first term generates the constituent quark mass and the second one gives rise to a one- 
boson exchange interaction between quarks. The main contribution of the third term comes 
from the two-pion exchange which will be simulated by means of the one-sigma exchange 
potential. Based on the nonrelativistic reduction of the former Lagrangian, one obtains an 
interaction between quarks mediated by the exchange of different Goldstone bosons. Among 
them, the one-pion, one-kaon and one-eta exchanges contribute with a central and a tensor 
interaction, while the one-sigma exchange gives a central and a spin-orbit term. Explicit 
expressions of these potentials can be found elsewhere [6] . The parameters of the Goldstone 
boson exchange interaction are fixed assuming that SU(3) flavor-symmetry is exact, only 
broken by the different mass of the strange quark. In the heavy quark sector, chiral symmetry 
is explicitly broken and therefore these interactions will not appear. 

For higher momentum transfer quarks still interact through gluon exchanges. Following 
de Rujula et al. [7] the one-gluon-exchange (OGE) interaction is taken as a standard color 
Fermi-Breit potential. In order to obtain a unified description of light, strange and heavy 
mesons a running strong coupling constant has to be used [4]. The perturbative expression 
for a s (Q 2 ) diverges when Q — > Aq C d and therefore the coupling constant has to be frozen at 
low energies. We parametrize this behavior by means of an effective scale dependent strong 
coupling constant [8] 
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where \i is the reduced mass of the qq system and ao, /io and A are fitted parameters [6]. 
This equation gives rise to a s ~ 0.54 for the light quark sector, a value consistent with 
the one used in the study of the nonstrange hadron phenomenology [9], and it also has an 
appropriate high Q 2 behavior, a s ~ 0.127 at the Z mass [10]. The S function appearing in 
the OGE has to be regularized in order to avoid an unbound spectrum from below. Taken 
into account that for a coulombic system the typical size scales with the reduced mass, we 
use a flavor-dependent regularization ro(/x) = ro/fj, [6]. Moreover the Schrodinger equation 
cannot be solved numerically for potentials containing 1/r 3 terms. This is why we regularize 
the noncentral terms of the OGE with a monopole form factor with parameter A g (/x) = /x/A° 

[6]- 

The other nonperturbative property of QCD is confinement. Lattice QCD studies show 
that qq systems are well reproduced at short distances by a linear potential that it is screened 
at large distances due to pair creation [11]. One important question which has not been 
properly answered is the covariance property of confinement. While the spin-orbit splittings 
in heavy quark systems suggest a scalar confining potential [12], in Ref. [13] showed that 
the Dirac structure of confinement is of vector nature in the heavy quark limit of QCD. On 
the other hand, a significant mixture of vector confinement has been used to explain the 
decay widths of P-wave D mesons [14]. Such property being irrelevant for the central part of 
the interaction, determines the sign and strength of the spin-orbit Thomas precession term. 
Therefore, we write the spin-orbit contribution of the confining interaction as an arbitrary 

combination of scalar and vector terms V^oNi^ij) = (1 — a s)Vv°(^ij) + a s^s°{^ij) where 
ySO^\ySO(^f^ j g ^ e vector (scalar) spin-orbit contribution. 

III. RESULTS 

With the quark-quark interaction described above we have solved the Schrodinger equa- 
tion for the different qq systems. Let us briefly discuss the parameters of our model. Most 
of the parameters of the Goldstone boson fields are taken from the NN sector. The eta 
and kaon cutoff masses are related with the sigma and pion one as explained in Ref. [15]: 
A[u(<i)s] ~ A(ud) + m s where m s is the strange quark current mass. The confinement pa- 
rameters a c and fi c are fitted to reproduce the energy difference between the p meson and its 
first radial excitation and the J/ip and the ip(2S). The parameters involved in the OGE are 
obtained from a global fit to the hyperfine splittings well established in the Particle Data 
Group (PDG) [16]. Finally, one has to fix the relative strength of the scalar and vector 
confinement. Using a pure scalar confining potential one obtains 1363 MeV and 1142 MeV 
for the ai(1260) and a 2 (1320), respectively, in complete disagreement with the order and 
magnitude of the experimental data. Introducing a small mixture of vector confinement, 
a s =0.777, the experimental order is recovered, being now the masses 1198 MeV and 1322 
MeV, respectively, both within the experimental error bars. This value also allows to obtain 
a good agreement with the experimental data in the cc and bb systems. 

In Fig. 1 we show results for the light pseudoscalar and vector mesons and for heavy 
mesons. The agreement with experimental data is remarkable. Let us emphasize that with 
only 11 parameters we are able to describe more than 110 states [6]. 
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A. The rj c (2S) and D SJ states 



Recently Belle and BaBar collaborations have reported new experimental measurements 
that immediately prompted different interpretations. The Belle collaboration has reported 
the following values for the mass of the rj c (2S) 



that are significantly larger than most predictions of constituent quark models and the 
previous experimental value of the PDG: M[r] c (2S)} = 3594 ± 5 MeV. It has been pointed 
out that these values cannot be easily explained in the framework of constituent quark models 
because the resulting 25' hyperfine splitting (HFS) would be smaller than the predicted for 
the IS ones. In fact the predicted ratio of the 2S to IS* HFS in charmonium is 



whereas the experimental one is R=0.273 if M[r] c (2S)} = 3654 MeV, R=0.547 for 3622 MeV, 
and R=0.479 for 3630 MeV. In view of these differences some authors have claimed for an 
a s coupling constant depending on the radial excitation. 

Our result is M[r) c (2S)] = 3627 MeV, within the error bar of the last two Belle measure- 
ments, the ones obtained with higher statistics. Moreover the ratio 2S to IS HFS is found 
to be 0.537, in perfect agreement with the experimental data. The reason for this agreement 
can be found in the shape of the confining potential that also influences the HFS, the linear 
confinement being not enough flexible to accommodate both excitations [23] . 

The experimental measurement reported by BaBar is a narrow state near 2317 MeV 
known as -0^(2317) [24]. This state has been confirmed by CLEO [25] together with 
another possible resonance around 2460 MeV. Both experiments interpret these resonances 
as J p = + and 1 + states. This discovery has triggered a series of articles [26] either 
supporting this interpretation or presenting alternative hypothesis. 

The most striking aspect of these two resonances is that their masses are much lower 
than expected. The mass difference between the PDG mass value _Di(2420) and the lowest 
value of Dq reported by the Belle collaboration is about 135 MeV. Using these data and the 
PDG mass value for the -0^(2536) one would expect the D Sq mass to be around 2400 MeV, 
almost 100 MeV greater than the measured values. Our results are shown in Table I. They 
agree reasonably well with the values of the PDG for both for the D's and D s 's states, but 
fail to reproduce the two D s states reported by Belle and CLEO. We also agree with the 
results by FOCUS [28] for the D state but we are far from the result of Belle for the same 
state. 



It is still not clear which are the members of the ++ nonet corresponding to L — S — 1 
qq multiplets. There are too many ++ mesons observed in the region below 2 GeV to be 



' 3654 ± 6 ± 8 MeV Ref. [17] 
M[r] c (2S)} = I 3622 ± 6 ± 6 MeV Ref. [18] 
( 3630 ± 8 MeV Ref. [19] 



R 



AM* FS 
AMfi FS 



' 0.84 Ref. [20] 
< 0.67 Ref. [21] 
0.60 Ref. [22] 



IV. THE SCALAR SECTOR: qq STUDY 
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explained as qq states. There have been reported in the PDG two isovectors IJ PC = 10 ++ : 
ao(980) and a (1450); five isoscalars IJ PC = 00++: / (600), / (980), / (1370), / (1500) and 
/o(1710); and three IJ PC = ±0++: X o *(1430), fT *(1940) and recently «(900). The naive 
quark model predicts the existence of one isovector, two isoscalars and two 1=1/2 states. 

Our results are shown in Table II. Using this table one can try to assign physical states 
to ++ nonet members. We observe that there are no nonstrange states with 1 = and 
mass close to 1 either 1.5 GeV, which would correspond to the /o(980) and the /o(1500). 
The same occurs in the strange sector around 0.8 GeV, we do not find a qq partner for the 
«(900). 

Let us discuss each state separately. With respect to the isovector states, there appears a 
candidate for the ao(980), the 3 Pq member of the lowest 3 Pj isovector multiplet. The other 
candidate, the ao(1450), is predicted to be the scalar member of a 3 Pj excited isovector 
multiplet. This reinforces the predictions of the naive quark model, where the LS force 
makes lighter the J = states with respect to the J = 2. The assignment of the a (1450) 
as the scalar member of the lowest 3 Pj multiplet would contradict this idea, because the 
a2(1312) is well established as a qq pair. The same behavior is evident in the cc and the bb 
spectra, making impossible to describe the ao(1450) as a member of the lowest 3 Pj isovector 
multiplet without spoiling the description of heavy-quark multiplets. However, in spite of the 
correct description of the mass of the a (980), the model predicts a pure light-quark content, 
what seems to contradict some experimental evidences. The a (1450) is predicted to be also 
a pure light quark structure obtaining a mass somewhat higher than the experiment. 

In the case of the isoscalar states, one finds a candidate for the /o(600) with a mass of 402 
MeV, in the lower limit of the experimental error bar and with a strangeness content around 
8%. The / (980) and / (1500) cannot be found for any combination of the parameters of the 
model. It seems that a different structure rather than a naive qq pair is needed to describe 
these states. The / (1500) is a clear candidate for the lightest glueball [29] and our results 
support this assumption. Concerning the /o(1370) (which may actually correspond to two 
different states [30]) we obtain two almost degenerate states around this energy, the lower 
one with a predominantly nonstrange content, and the other with a high ss content. Finally 
a state corresponding to the / (1710) is obtained. 

Concerning the 1 = 1/2 sector, as a consequence of the larger mass of the strange quark 
as compared to the light ones, our model always predicts a mass for the lowest ++ state 
200 MeV greater than the ao(980) mass. Therefore, being the ao(980) the member of the 
lowest isovector scalar multiplet, the k(900) cannot be explained as a qq pair. We find a 
candidate for the 1^(1430) although with a smaller mass. 

In conclusion our results indicate that the light scalar sector cannot be described in a 
pure qq scheme and more complicated structures or mixing with multiquark states seems 
to be needed. Our conclusion concerning the /o(980) and the k(900) are very similar to 
the ones obtained by Umekawa [31] using the extended Nambu-Jona-Lasinio model in an 
improved ladder approximation of the Bethe-Salpeter equation. This fact seems to indicate 
that relativistic corrections would not improve the situation and the conclusions remain 
model independent. 



5 



V. THE SCALAR SECTOR: TETRAQUARK STUDY 



In the naive quark model, to construct a positive parity state requires a unit of angular 
momentum in a qq pair. Apparently, this takes an energy around 1 GeV since similar meson 
states (1 ++ and 2 ++ ) lie above 1.2 GeV. However a more complicated structure, like q 2 q 2 , 
suggested twenty years ago by Jaffe [32] can couple to ++ without orbital excitation and 
therefore could be a serious candidate to explain the structure of some light scalar mesons. 

In this section we study tetraquark bound states, focusing our attention in those states 
with the quantum numbers of the scalar mesons. We solve the Schrodinger equation using a 
variational method where the spatial trial wave function is a linear combination of gaussians 
The technical details are given in Ref. [33]. Let us only mention that neglecting the exchange 
terms in the variational wave function is fully justified for the heavy-light tetraquarks but it 
may induce some corrections in the present study. Due to the presence of the kaon-exchange 
there is a mixture among different configurations with the same isospin. In particular, in 
the isoscalar sector the configurations: [(qq)(qq)], [(qs)(qs)], and [(ss)(ss)] are mixed. The 
same happens in the isovector case for the configurations: [(qq)(qq)], and [(qs)(qs)], and in 
the I = 1/2 case for the configurations: [(qq)(qs)], and [(gs)(ss)]. In all cases q stands for a 
u or d quark. 

The results are shown in Fig. 2. We present the lowest states for the three isospin sectors. 
The x's show the results obtained within the model described in Sect. II. As one can see, 
there appear two states, in the isoscalar and isovector sectors, with almost the same mass, 
although too high to be identified with the /o(980) and ao(980). In the 1=1/2 sector, 
there appears a candidate to be identified with the k(900). It has been recently argued 
the possible importance of three-body forces arising from the confining interaction for those 
systems containing at least three quarks [34]. We have performed a calculation including a 
three-body confining term as the one reported in Ref. [34]. The strength of this interaction 
has been chosen to reproduce the mass of the /o(980). The results are plotted in Fig. 2 by 
*'s. As one can see the degeneracy between the isoscalar and isovector states remains, while 
the lowest state of the isoscalar and 7=1/2 sectors are almost not affected. From these 
results one could inferred that the scalar sector needs the presence of tetraquark structures 
to be understood. 

As a summary, we have obtained candidates for all the light scalar mesons reported in 
the literature. Our results suggest that there would be some states, as it is the case of 
the ao (980) and /o(600) that would be a mixture of a qq and tetraquark structure, but it 
definitively assigns a tetraquark structure to the / (980) and the k(900). 
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TABLES 



TABLE I. D's and D s 7 s masses in MeV compared to all known experimental data [27] 





M*(0+) 


M(l+) M(l+) 


M*(2+) 


M(0")(2S) 


M*(1-)(2S) 


£>'s 


PDG (0) 




2422±2 


2459±2 






PDG (±) 




2427±5 


2459±4 


2637±7 


2637±7 


FOCUS (0) 


~ 2420 




2463±2 






FOCUS (±) 


~ 2420 




2468±2 






Belle (0) 


2290±30 


2400±36 2424±2 


2461±4 






CLEO (0) 




2461±51 








Our result 


2437.9 


2496.0 2495.7 


2500.5 


2641.7 


2700.1 


D s 's 


PDG (±) 




2535.3±0.6 


2572.4L1.5 






FOCUS (±) 




2535.1±0.3 


2567.3L1.4 






BaBar 


2317 










CLEO 


2317 


2460 








Our result 


2470.5 


2565.5 2549.7 


2584.4 


2698.7 


2764.3 


TABLE II. Light scalar meson masses in MeV 


State (n 2I+1 > 2S+1 Lj) 


Meson 


Our result 


Experiment 


1 3 ' 3 P 




a (980) 


983.5 




984.7L1.2 


2 3 > 3 P 




a (1450) 


1586.3 




1474±19 


l^Po 




/o(600) 


402.7 




400-1200 


l^Po 




/o(1370) 


1341.7 




1200-1500 


2 1 ' 3 P 




/o(1370) 


1391.2 




1200-1500 


2 1 > 3 P 




/o(1710) 


1751.8 




1713L6 


3 1 ' 3 P 




/o(2020) 


1893.8 




1992±16 


3 1 ' 3 P 




/o(2200) 


2212.2 




2197L11 


l 2 ' 3 Po 




#5(1430) 


1213.5 




1412L6 


2 2 ' 3 P 




# *(1950) 


1768.5 




1945±30 
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FIG. 1. Meson spectra (a) light peudoscalar and vector mesons (b) bb mesons 
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FIG. 2. Tetraquark masses for the scalar sector 
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